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Q Abstract 

We formulate a Euclidean lattice theory of interacting elementary spin-half elec- 
tric and magnetic charges, which we refer to as electrons and magnetic monopoles 
respectively The model uses the polymer representation of the fermion determi- 
nant, and exhibits a self-dual symmetry provided electric charge e and magnetic 
charge g obey the minimal Dirac quantisation condition eg = 2n. In a hopping 
parameter expansion at lowest order, we show that virtual electron and monopole 
loops contribute radiative corrections of opposite sign to the photon propagator. 
We argue that in the limit e — > 0, fermion mass \x — > 0, the model describes QED 
together with strongly interacting monopoles whose chiral symmetry is sponta- 
neously broken. Prospects for the existence of an interacting continuum limit at 
the self-dual point e = g are discussed. 
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1. Introduction 

Ever since Maxwell's equations were written down, it has been natural to specu- 
late why there exist fundamental electric charges but apparently no magnetically-charged 
counterparts. At the classical level, it is straightforward to modify Maxwell's equations for 
electric and magnetic field strengths (in appropriate units) by the introduction of magnetic 
charge and current densities (p m ,j m ) analogous to the electric quantities (p e , j e ): 

<9 E = VAB-j e , 
8 B = -VAE-j m , 
V.E = Pe , 
V.B = p m . 

The system of equations (1.1) is invariant under a duality transformation: 

E — E ; B ^ -E 

(1.2) 

(PeJe) >-»• (PmJm) 5 (Pmjm) ^ "(Peje)- 

The cost of this modification is that it is no longer possible to describe the electromagnetic 
field in terms of an abelian vector potential which is both globally defined and non- 
singular. This raises difficulties, because the standard quantum-mechanical description of 
electromagnetism is formulated in terms of the potential rather than the field strengths, 
and indeed the success of the prototype quantum field theory, QED, depends critically on 
the local gauge symmetry transformations which relate physically equivalent configurations 
of the field. Dirac [1] showed how to introduce magnetic monopoles as pointlike sources 
of magnetic flux by connecting them to line singularities in A^; so-called Dirac strings. 
Despite the singularity the string has no observable effect on a particle carrying electric 
charge e, in either classical or quantum mechanics, if the strength of the magnetic charge 
g is quantized according to the famous condition 

eg = 2mrh, (1-3) 

where n is an integer. Planck's constant will henceforth be set to one. The spatial position 
of the string can be changed under gauge transformations, but the magnetic flux emanating 
from the monopole is uneffected. The Dirac quantization condition (1.3), or variations of 
it, has been central to all subsequent attempts to create a quantum theory which includes 
monopoles (for a comprehensive review see [2]). 
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Attempts to construct a quantum field theory containing both particles with electric 
charge and particles with magnetic charge - so-called Quantum ElectroMagnetoDynamics 
(QEMD) - began with Cabibbo and Ferrari [3] and Schwinger [4]. The latter formulation 
is in the Hamiltonian framework, and contains fermionic electron and monopole fields each 
interacting with its own vector potential field. The two potentials are not independent. 
The model also has an explicit dependence on a vector s defining the direction of a Dirac- 
type singularity, needed to relate the potentials to the current distributions defined by the 
matter fields. Later an equivalent Lagrangian formulation was given by Zwanziger [5]. In 
both cases the quantization condition was modified to be 

eg = 4mr. (1.4) 

The string-dependence obscures the Lorentz invariance of the results. In addition, since 
the string's position is gauge-dependent, insistence on a fixed string vector s M leaves no 
possibility of a gauge invariant formalism. An alternative approach, beginning with [3], 
allows the string, rather than being fixed and infinite, instead to run between monopole 
- antimonopole pairs, resulting in a description in terms of closed particle world-lines 
spanned by surfaces known as Dirac sheets. This formulation is gauge invariant and stresses 
the path-dependence of the particle dynamics - and hence is inherently non-local. The 
requirement of surface-independence once again results in the Dirac quantization condition. 
The world-line approach was used by Brandt, Neri and Zwanziger [6] to demonstrate that 
gauge-invariant Green functions in the Lagrangian formulation are independent of the 
vector defining the string. 

As we have sketched, a common feature of the various formulations, some of which 
have been shown to be equivalent [2], has been the impossibility of constructing a man- 
ifestly local, covariant and gauge invariant formulation. Another problematic issue has 
been how to proceed once the model is formulated - various authors have succeeded in 
finding "Feynman rules", but these are of limited use due to the condition (1.3): QEMD 
is strongly interacting. Indeed, the origin of the quantization condition, the requirement 
that exp(ieg) = 1, is inherently non-perturbative. 

Calucci and Jengo [7], realising that the strongly-interacting nature of the problem 
necessitates new techniques, regularised a version of the worldline formalism [8] on a Eu- 
clidean spacetime lattice. A non-perturbative regularisation must be used if the ultraviolet 
behaviour of the model is to be discussed adequately. Their model has two field strength 
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tensors, defined on the direct and dual lattices respectively, and requires the more restric- 
tive quantisation condition (1.4). The lattice is used as little more than a formal framework, 
but enables the use of a proper time technique to calculate the effects of virtual particle 
loops on the interaction among electric and magnetic charge strengths. Their claim is that 
if eo and go are "bare" charge strengths satisfying (1.4), then the "renormalised" charges 
cr, 9r satisfy 



ie, eg is a renormalisation group invariant. The result follows essentially, as we shall 
see, because of an extra factor of i in the electron - monopole interaction which reverses 
the sign of the vacuum polarisation diagram: monopoles "anti-screen" electric charge, 
and vice-versa. Although persuasive, this result is in contradiction to that obtained by 
other authors [9]. Moreover, the Feynman rules developed by Calucci and Jengo contain an 
ordering ambiguity which requires an artificial distinction to be drawn between valence and 
virtual current loops, once again reflecting the non-perturbative nature of the quantisation 
condition and suggesting that any approach based on series expansions in e and g will not 
be fully satisfactory. 

As pointed out in [7], the self-dual point e = g appears to define a fixed point of 
the renormalisation group, since if the duality symmetry is respected by the regularisation 
then en = gn and engn = constant together imply that neither e nor g is renormalised. 
The existence of an ultraviolet fixed point, or continuum limit, for abelian theories such 
as QEMD or indeed QED has been a matter of speculation for many years [10]. In 
perturbative QED we know the relation between bare and renormalised electric charge: 



where A is the cutoff and fx a physical scale, implying that — > as A/ fx — > oo. This 
phenomenon is known as triviality - QED can only exist as an interacting theory if the 
cutoff is finite. An alternative scenario requires the existence of a zero in the /3-function 
describing renormalisation group flow of the interaction strength beyond the perturbative 
regime. One possibility was advanced by Miranskii and co-workers, who investigated the 
self-energy of the electron using a truncated system of Schwinger-Dyson equations [11]. For 
a critical value e c = 2^/^/3 the chiral symmetry protecting the electron from acquiring 
mass is spontaneously broken. A physically meaningful ground state can only be found if 



dRgR = e g , 



(1.5) 



eg(A) 



(1.6) 
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the bare charge e obeys 



eg (A) 



) 



2 



= 1 + 



7V 



(1.7) 



ln(A//i) 



where /U is now related to the electron mass in the broken phase. Hence the continuum 
limit is taken as e — > e c . Note that e c is greater than either of the "minimal" self-dual 
values in QEMD e = V2tt, y/An. 

The Miranskii fixed-point hypothesis received support from numerical simulations of 
non-compact lattice QED (NCQED) [12], where it was found that for sufficiently strong 
electric charge the model exhibits a continuous phase transition from the perturbative 
phase to one in which chiral symmetry is broken, signalled by the non-vanishing of the 
condensate {ipif})- However, intensive computer simulation has not found the scaling be- 
haviour consistent with Miranskii's prediction of an essential singularity as e — » e c +: al- 
though there has been some dispute in the literature, currently the existing numerical data 
for the equation of state relating (V>V)> e an d the bare electron mass \i are best fitted by a 
critical coupling e c in the region 2.2-2.3 (for the case of four dynamical fermion species), 
that is, less than the value predicted by the Scwinger-Dyson approach, with a power-law 
singularity with either non-mean field exponents [13], or mean field exponents with log- 
arithmic corrections [14]. In either case the the continuum limit is thought to consist of 
both electrons and bosonic e + e~ bound states. The mean field fit is motivated by a picture 
in which the bound state charge radii and physical interaction strengths vanish so that the 
model is trivial. 

To a large extent, the conclusions drawn from numerical simulations are subject to 
theoretical prejudice; the true fixed point theory, if it exists, may contain many other inter- 
actions apart from the electron-photon interaction of the minimal gauge coupling (alterna- 
tives for the lattice electron-photon interaction have been explored in [15,16]). Constraints 
on interaction Lagrangians due to the requirement of perturbative renormalisability may 
not apply at a strongly-coupled fixed point. Operators such as the chirally-symmetric 
(ifji/;) 2 — (■075'0) 2 may acquire large anomalous scaling dimensions through non-perturbative 
effects and hence become relevant at the fixed point [17]. It is difficult to know a priori 
how to restrict the space of bare theories which need to be explored before an interacting 
continuum limit can be excluded. A defensible viewpoint is that the mere existence of an 
observed regime where perturbative QED applies (as an effective theory) implies the exis- 
tence of an underlying theory containing massive degrees of freedom and/or short-ranged 
interactions which are so far imperfectly determined. Of course, it is more usual to consider 
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electromagnetism as part of a unified description of fundamental interactions, though it is 
worth recalling that the Standard Model gauge group contains U(l) as a subgroup, so the 
issue of a continuum limit for abelian gauge theories cannot as yet be ignored. 

In this paper we shall explore the notion that the bare theory, and hence the fixed- 
point theory, are constrained to have the symmetry (1.2) of electromagnetic duality, by 
modifying lattice QED so that it has both electron and fermionic magnetic monopole 
degrees of freedom; in other words, we construct a lattice QEMD. We have been motivated 
in part by the observation that at the large values of e in the vicinity of the observed chiral 
transition, the lattice photon field may be interpreted as containing line dislocations defined 
on the dual lattice, which are akin to magnetic currents [18]. The monopole dynamics 
in NCQED are unconventional, however, since the Dirac sheet in this formulation has 
non-vanishing action density, which means that in effect the monopoles probably exist in 
tightly bound dipole pairs with no associated long-ranged electromagnetic fields. Lattice 
monopoles bear a closer resemblance to Dirac ones when they occur as excitations in 
compact or periodic abelian gauge theory [19,20], in which the Dirac sheet costs no action 
and is physically unobservable; in this case the monopoles act as sources for long-ranged 
Coulomb fields. Banks et al [19] showed how the partition function for the Villain form of 
U(l) lattice gauge theory may be represented as a Coulomb gas of monopole point charges 
in three Euclidean dimensions or monopole world line loops in four dimensions by using an 
exact lattice duality transformation. By combining this formulation of lattice monopoles 
with the ideas of Calucci and Jengo, we are able to construct a theory with electron and 
monopole currents interacting via a Coulombic potential, with the charge strengths given 
by the less restrictive quantisation condition (1.3) with n = 1. Field strengths and gauge 
potentials, with all their associated complications [2], play a secondary role. The cost is 
that, as in all previous attempts, we are unable to give a local covariant action. Instead, 
there is a partition function written in terms of "polymer configurations" on both direct 
and dual lattices, based on the polymer representation of the lattice fermion determinant 
introduced by Karowski et al [21]. Because the fundamental excitations are fermions, an 
additional penalty is that the terms in the partition function fluctuate in sign. In the limit 
e — > 0, g — > oo, we argue that the electron and magnetic sectors decouple, and that the 
model reduces to a local action describing electrons interacting via photon exchange as 
in QED. The monopoles' chiral symmetry is spontaneously broken, so that the only light 
degrees of freedom in the monopole sector are monopole - antimonopole "magnetopions" 
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corresponding to the Goldstone mode associated with the broken symmetry. The only 
conceivable point where both electrons and monopoles can exist as light degrees of freedom 
is the self-dual point, which is therefore a candidate continuum limit for the model. 

The rest of the paper is organised as follows. In Section 2 we introduce lattice QEMD 
in a heuristic fashion, by first discussing non-compact lattice QED in the polymer represen- 
tation, and then reviewing the transformation from Villain lattice model to the monopole 
gas, before finally writing the QEMD partition function. We also briefly review connec- 
tions with other approaches to QEMD. In section 3 we discuss the possible phase structure 
of the model, using exact results, analytic approximations and known numerical results. 
A hopping parameter expansion to lowest order is used to find that electron and monopole 
loops induce radiative corrections of opposite sign in the photon propagator, in agree- 
ment with [7]. Unfortunately it seems that in the minimal model there is no value of the 
couplings, including the self-dual point, for which both electrons and monopoles are light 
degrees of freedom; therefore probably no interacting continuum limit exists. The reasons 
are discussed in the final section, along with possible modifications of the model which 
may improve the prospects for a continuum limit, and phenomeno logical consequences. 

2. A Lattice QEMD 

In this section we build the QEMD partition function in simple stages. We begin with 
the partition function of non-compact lattice QED (NCQED): 



where the determinant results from integrating over fermion fields of bare mass fi, leaving 
the photon fields 9^ G (— oo, +oo) as real dynamical variables defined on the links of a 
hypercubic lattice. The field strength tensor Mi , is then defined by 



with A+ the forward lattice difference operator. We choose to use the staggered fermion 
formulation (this is not crucial in what follows), so the fermion kinetic operator is given 




(2.1) 



(2.2) 



by 



P( x > V) = g XI [<W+A exp(ie9^(x)) - 5 y)X _ A exp(-ie9^(x -£))], (2.3) 
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where the Kawamoto- Smit phases rj^x) = (— l) XlH i~ x m-i. The electric charge e enters 
only via the gauge covariant connection in the derivative. To avoid a proliferation of sum- 
mation signs in the formulas which follow we adopt the convention that repeated spacetime 
indices and repeated spatial arguments are summed over, even in squared quantities. We 
shall also employ the subscript £ in connection with oriented variables associated with a 
particular link. 

The partition function in (2.1) is based on an action with a local gauge invariance, of 
course. For Z to be strictly well-defined the gauge action must be gauge fixed to avoid the 
partition function from diverging on integration over a flat direction. We now immediately 
introduce one of the major conceptual tools we will need - the polymer representation of 
the fermion determinant introduced by Karowski et al [21]. We write 

Z = J V9» N |u^^phase(C^)|exp [-^(x)) , (2.4) 

where k is the hopping parameter 1/2// and C is a partition of the lattice in which every 
site is either isolated (a monomer), joined to a nearest neighbour by a double "bond" (a 
dimer) , or visited once and once only by a closed oriented self-avoiding loop of single bonds 
(a polymer). The set of polymers in a particular configuration C will be denoted {T}. Each 
configuration C corresponds to a unique term in the expansion of the determinant: the N m 
monomers correspond to diagonal elements, and hence contribute a factor fjL Nm , and the 
Ng bonds correspond to off-diagonal hopping terms. Clearly we have 

N m + N e = N (2.5) 

where N is the lattice volume. The phase factor has a number of different components: 
a signed permutation factor from the expansion of the determinant; a sign from the path 
over forward and backward links and the Kawamoto-Smit phases (we shall refer to this 
as the "signed KS phase"); and finally the gauge connection n^ec exp(ie^). A moment's 
thought reveals that monomers and dimers give trivial contributions to the phase - only the 
polymers {V} contribute. The smallest polymer which gives a negative sign is a non-planar 
loop of length six. We can rewrite Z as follows: 



Zoc fveY, k~ n -(-1) n - f II A I II U{T)\ew(-\*U*) 




re{r} 

Zphoton P] j 



(2.6) 



where N r is the number of polymers, r\g is the signed KS factor associated with the link t, 
and U(T) = n^er ex P(^ e ^) ^ s tne Wilson loop associated with the polymer T (note that 
since the polymers are oriented, then both U(T) and U*(T) appear as separate terms in Z 
with equal weight). Z p hoton[F] is simply the partition function for free photons coupled to 
a Wilson loop distribution {]?}. Note the sign factor associated with the number of loops: 
the relation with Fermi-Dirac statistics is clear. Note also that despite the superficial 
similarity this is not the conventional hopping parameter expansion - the sum over {C} is 
finite on a finite system, and this expression for Z is valid for all \i. 
Now let's focus on Z p h oton \T\. This can be written 

Z photon [T]= f + UdOeexpl-jiA+OvW-A+O^f + ie ]T j e e \ , (2.7) 
J -°° e \ ee{r} J 

where jg is the characteristic function of the polymer configuration, ie 

{+1, if £ a forward link e {T}; 
-1, if £ a backward link e {r}; (2.8) 
0, otherwise. 

Of course, jg is nothing other than the electric current. By performing a shift of the 
integration variables 0g, we can rewrite 

z P hoton[j] = z photon [o] exp ^-yj^Wvl^ - y)jv{y) s j ■ (2.9) 

Here v^ u (x) is the lattice Coulomb propagator, whose existence requires a gauge-fixing 
term in the photon action. In Feynman gauge it satisfies 

A+A-v ltu (x) = -6 liV 6(x), (2.10) 

where A~ is the backward difference operator. We have written the arguments x, y to 
stress that the interaction between the current loops is non-local in this representation. The 
lattice Coulomb propagator is finite for zero spatial separation, and for large \x\ behaves 
as |x| -2 . Hence, in the polymer representation of the fermion determinant, lattice QED 
resembles a Coulomb gas of unit charged non-intersecting loops of variable sign: 



Zncqed oc Z photon [0]J2^~ Nm (-^) Nr I II ^ I ex P ( J2 -\wu>3t' ] • (2.11) 

{C} ye{r> / V,fe{r} 



It is important to note, however, that the role of monomers and dimers is also important, 
as we shall see. 

Next, we review another well-known model - the Villain approximation to U(l) lattice 
gauge theory [19]. 

n ^ exp \ 2 n,tu ^ Q ' iu ^ ~ ~4 n l^) ' ( 2 - 12 ) 

where n^ v {x) are integers defined on each plaquette of the lattice. The Villain action 
looks a little strange, but is simply the convolution of the Gaussian non-compact gauge 
action 2 with an infinite array of delta functions having period 2tv. Around any one 
of these locations the Villain action has approximately the Gaussian form of NCQED, 
but by construction it is also periodic in 6, and in this respect is similar to the Wilson 
form (1 — cos©). Its great virtue is the existence of a sequence of exact transformations 
which mean its phase structure may be understood in terms of monopole excitations. The 
integral over {9} can be performed immediately to leave a constrained action of the n^ u : 

Zvillain = ^ <H A *>^) eX P (-^T n l,s( X )) • ( 2 - 13 ) 

w v J 

The constraint can be solved by rewriting in terms of another integer variable /: 

n^(x) = e^^Ap^x). (2.14) 

Here, / is defined on the links of the dual lattice, as indicated by the tilde on x, and e is 
the totally antisymmetric tensor. As a rule, any expression containing an odd number of 
e symbols relates variables living on direct and dual lattices. Note that (2.14) does not 
uniquely specify /, or even require it to be integer- valued, since the definition of n remains 
invariant under 

l^x)^l^x)+A+A(x), (2.15) 
where A is an arbitrary scalar function. Now we can write 

Zvillain = $>xp (~(A+l v {x) - A+l^x)) 2 ) . (2.16) 

{1} v J 

This expression is very reminiscent of NCQED, except that the dynamical variables are 
integers, not real numbers. The partition function may now be reexpressed, first by using 
the Poisson formula 

/oo 
#/(</>)e 2 ™^, (2.17) 
-oo 



{1} {m} 



10 



where the m are integer, and then by shifting the variable in the resulting Gaussian 
integral as before. The result is 



Zvillain = Zphoton ^ eXP ( S ~Y m t°U' m i' I • ( 2 - 18 ) 
{m} \W ~ J 

Here, are integer variables defined on the links of the dual lattice, which obey the 
constraint A~m^(x) = (since otherwise (2.18) would not be invariant under gauge 
transformations (2.15) on 0). These two facts constrain the m to form closed loops on 
the dual lattice, which interact via the Coulomb potential as in NCQED. Once again the 
partition function may be thought of as a sum over polymer configurations, though this 
time the polymers need be neither singly charged nor self-avoiding. The coupling strength 
is g = 2%/e; the excitations thus resemble magnetic current loops with charge specified by 
the quantisation condition (1.3). 

We can confirm this picture of interacting magnetic current loops, by coupling gauge- 
covariant fermions to the Villain model, to yield a model we will call "periodic QED" 
(PQED) to distinguish it from NCQED. From what we already know, we can write 

Zpqed = J V9 det(p[9] + fi)^exp Qn^(x)0^(x) - e —n\ v {x) 



{n} 




(2.19) 



Here, Zvui a in[^} is the same partition function as that of (2.12), but this time in the 
presence of an electric current distribution defined by {r}. The steps (2.12-18) can be 
repeated to yield 



x 

{m} 



zp QED =z photon [o]j2*- Nm (-v Nr \ n w EtaEUs. 

{C} Ve{r} J {j} {S} 

exp ( - yJ/*( a O u /«'( a; ~ y)jv(y)^J ex P y"v(£)*v( £ ^ y) m v(y)^J 

x exp ^-2Trim^(x)v^(x - y)^e I/Ap(J A+5' p(J (y + u^j . 

(2.20) 

To implement the steps (2.12-18) it has been necessary to introduce the oriented char- 
acteristic function ,S Mi , for a surface which spans the self-avoiding loops in the polymer 

11 



expansion, but is otherwise unconstrained. Due to the Kronecker <5's both j and S may 



However, in the form (2.20) the similarity between the electric and magnetic currents is 
suggestive. The model has a new term - the last exponential in (2.20) - which describes 
the interaction between electric and magnetic currents (the spacetime argument of S is 
determined by the requirement that the Coulomb propagator couples to a conserved cur- 
rent at either end). To see that this is indeed an interaction between magnetic and electric 
current distributions, consider the distribution of a static magnetic charge mo(x) with a 
spacelike surface Sij(y + 0), with i,j ranging from 1 to 3. The interaction may be written 
in the form igmo(x)V mag (x) where V mag is the magnetic scalar potential: 



The factor of % is a consequence of the formulation in Euclidean space. In the long wave- 
length limit the contribution to Vm ag from a current loop spanned by a surface S assumes 
its textbook form: 



where O is the solid angle subtended by the circuit at x. Hence the interaction only depends 
on the boundary of the loop dS. The magnetic scalar potential is not single valued, since 
O is only defined modulo 4%: however, this has no physical consequence in the present 
case due to the Dirac quantisation condition, as we shall now discuss. 

To see explicitly that the interaction only depends on the edge of the surface, consider 
the following sequence of transformations, called the "disentangling theorem" in ref. [8], 
on the monopole-electric current interaction term in Feynman gauge; we define to be 
the integer-valued characteristic function for a Dirac sheet which spans the m loop as S 



t This is not true in a finite volume - we are grateful to John Stack for pointing this 
out 



be regarded as auxiliary fields: a given configuration is uniquely specified by C and m 1 . 




k&t v 3d(x - y)S jk (y + 0) 



(2.21) 




(2.22) 
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spans the j loop: 

m^{x)v^{x - y)^e uXpa A+S pa (y + v) 

I 

= A~D Tp (x)v(x - y)-e flXpa A+S pa (y + p) 

= -^D Tp (x)A~A+v(x -y + r + p)e pXpcT S pa (y) 

= ^D Tp (x)A~A+v(x -y + f + p)e TpaX S p(7 (y) - ^D Tp (x)e Tppr7 S pr7 (y)5(x - y + r + p) 
= —^A~^D Tp (x -X-p- f)v(x - y)e TpaX A~ S pa (y) - ^D Tp (x)e Tpp(J S pr7 (x + t + p) 
= -Ja{y)v{y - x)]^e a \ Tp A+D Tp (x — X — p — r) — ^D Tp {x)e Tppa S pa {x + f + p). 

(2.23) 

We have used translation invariance v(x) = v(x) = v(—x), and the identity 

[e XTp aA-A+ + e pXp *A;A+ + e prXa A;A+ + e pTpX A- A+]v(x) = -e pTpa 5(x). (2.24) 

We have also made the lattice coordinates explicit, resulting in a slightly unwieldy expres- 
sion. This is due to the offset of the dual lattice origin from the direct one. More compact 
notations, which make use of the language of differential forms, are available (eg. [7,22]). 

The first thing to notice is that the expression \DeS is an integer (since D and S 
are themselves integer) - the intersection number of the two surfaces multiplied by the 
monopole and electron charges. Since the coefficient of the term is — 2ni = — ieg, this 
term contributes unity to Zpqed and hence has no dynamical influence. Therefore the 
electron-monopole interaction is manifestly independent of the interior of either of the 
surfaces S or D - which are unconstrained by the dynamics - but only their boundaries, 
respectively the electric current j or the magnetic current m. It is important to note 
that these currents are constrained to lie on direct and dual lattices, and that the surface 
independence, which is a necessary requirement for the recovery of a local limit for the 
quantum field theory, depends on the Dirac quantisation condition (1.3). Essentially the 
same argument was put forward originally in [6]. 

The interaction betwen electric and magnetic charges is best described on the lattice 
by a non-local interaction between conserved currents. However, it is useful at this stage to 
pause and review the connection between the formulation appearing in (2.20) and the more 
familiar description in terms of local potentials and field strengths. From the discussion 
preceding (2.21), we recall that magnetic charge has an interaction energy depending on 
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a magnetic scalar potential V mag . By treating V mag as the zeroth component of a four 
vector A^, it is possible to write the interaction term for electric current as ej^A^, and 
the interaction for magnetic current as gm^A^, where the two potential fields A^ and A^ 
are given, in continuum notation and in Feynman gauge, by 



-1 . .1 -1 



^ CfiXpcrOX Q2 '-'per- 



(2.25) 



Once again, the factors of % associated with the e symbols are due to the Euclidean space 
formulation. Field strength tensors F^ v and F^ v can now be defined in the usual way by 
taking the four-dimensional curl of the potentials; using (2.24) we find the relation between 
them: 

i ~ 
--e Mi ,A K (^A K + eS\ R ) = F^ v - gD^ v . (2.26) 

Hence both F and F contain ambiguities due to Dirac sheets, which in the continuum 
theory must be resolved by the introduction of compensation terms supported only on the 
sheets [2] . Of course, both surfaces S and D may be moved around at will by gauge trans- 
formations, and hence have no physical meaning. In continuum formulations of QEMD, in 
order to avoid double-counting, it is customary to express the ambiguity solely in terms of 
D^ v , but relations (2.25-26) expose the duality symmetry to the full. On the lattice prob- 
lems never arise so long as an action which is a periodic function of the field strength, such 
as (2.12), is chosen; the ambiguity then simply relates physically equivalent configurations 
of the F field (for the case of the Wilson action for U(l) lattice gauge theory see [20]). 

The final step in the construction of QEMD is to endow the monopole excitations 
with specific dynamical properties. In eqn. (2.20) the only difference between monopoles 
and electrons are the constraints imposed on the electric current loops by the polymer 
representation of the fermion determinant, ie. that C be a partition of the lattice into 
monomers, dimers and self-avoiding (and thus singly-charged) polymers. Suppose we now 
insist that the dual lattice be similarly partitioned by specifying C consisting of dual 
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monomers, dimers and polymers, which define loops of magnetic current. We then write 

z QEM D=z photon [o]J2Y K ~ Nm ~ N ' n (- 1 ) Nr+Nt ( II vt ) ( n m 

{C} {C} ye{r} J \Ie{f} 

X 5^^e{r} Y S j„,A^S^ Y 6 rn,le{r}Y 6 rn tl ,A^D^ 
{3} {S} {m} {D} 

x exp ("y^K^O* ~~ y^M^j ex P (-y m M( £ )*v( £ ^ y) m v(y)^ 

x exp ^27rzm M (5)^(5 - y^e^o-A^S^y + z>)^ . 

(2.27) 

Each configuration is completely specified by C, C; the variables j, m, S and D are auxiliary. 
Note that the monopole fields have been given the same hopping parameter k as the original 
fermion fields. Now, the crucial point is that the model defined by (2.27) is invariant under 
the following duality transformation: 

3 ix <-> nip e <-> g D ^ S Z <-> (2.28) 

The manipulations (2.23) play a crucial role in demonstrating this. Note further that 
under the interchange m M <-> — m M , every term in Zqemd remains unchanged (since the 
number of links in every polymer is even, the factors Y\l Ve remain unaltered) except the 
argument of the last exponential, which changes sign. Since this last term is a phase, and 
every m-polymer loop appears with equal weight going in either direction, we find that 

Zqemd = Zq EMD . (2.29) 

Therefore the model is self-dual, respecting the symmetry (1.2): in particular at the self- 
dual point e = g = V2n the electron and monopole degrees of freedom behave identically. 

To complete the formal exposition of the model, three final comments are needed. 
Firstly, although the role of the gauge potential has been suppressed in our treatment, 
notice that the interactions in the model (2.27) are all described in terms of conserved 
integer- valued currents; hence the model is clearly also gauge invariant. Next, although in 
the forms (2.6), (2.20) and (2.27) the various partition functions Zncqed, Zpqed and 
Zqemd consist of terms with fluctuating sign, in the first two cases of NCQED and PQED 
we know that the partition functions can also be expressed via a positive definite (albeit 
non-local) effective action, and hence are themselves positive definite on a finite lattice. 
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No such argument can be used for QEMD, since it is not expressible simultaneously in 
terms of both local electron and monopole fields. A demonstration of the positivity of 
Zqemd remains an open question, and is desirable for the discussion of the possible phase 
structure of the model (see next section). Thirdly, in the form (2.27) the magnetic current 
m M transforms as a vector, by analogy with the electric current. Even in the classical theory 
descibed by (1.1), it is impossible to satisfy invariance simultaneously under discrete parity 
(P) and charge-conjugation (C) transformations. For a vector m M , under P: 

Ek-E ; BhB; 

(2.30) 

(PeJe) !-»■ (fie, ~je) 5 (PmJm) ^ (Pm, ~jm), 

whereas under C: 

Eh-E ; B i-> -B; 

( 2 - 31 ) 

(Pe)je) l— ^ ( — Pei — je) j (Pmjjm) l— ^ ( — Pm, ~ jm)- 

The QEMD Maxwell equations (1.1) are invariant under C but not P, or hence CP. If 
is defined to be an axial vector, then under P: 

(PmJm) i-> (-PmJm), (2-32) 

and under C: 

), (2.33) 

in which case (1.1) are now invariant under P, but not C or CP. A Euclidean lattice 
analogue of this latter case can also be constructed if the monopole conserved current is 
derived using a kinetic term of the form 

(p + fi)(x,y) = \ y^r]^(x)e(x) [<y yia . +/i exp(ie^(x)) - 8 v , x -n exp(-ie# M (x - //)] 

2 m (2-34) 

+ fi6y jX e(x). 

Here, e(x) is the alternating phase (— i^i+^+a^+a;^ The result of this modification in the 
polymer expansion is to reverse the sign of every loop T containing (2n + 2) links. 



3. Phase Structure 

In this section we will try to discuss the lattice model described in section 2 more 
quantitatively, with particular emphasis on its phase structure. We shall draw insight 
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and information from what is already known about NCQED, Villain QED, and PQED. 
First, however, we shall consider the limit of large fermion bare mass fi, in which case an 
expansion in powers of 1/ ji is viable [23] . 

For large fx, the partition function in (2.27) is dominated by configurations in which 
N m and Nfh, the numbers of direct and dual monomers respectively, are large. Both dimers 
and polymers must be considered as excitations suppressed by powers of k = 1/2 fx. Only 
polymers, however, interact with the electromagnetic field, and hence with other loops, via 
the exponential terms of (2.27). Consider the smallest monopole loop excitation, simply 
four dual links bounding a dual plaquette. By considering the photon exchange diagrams 
of figure 1, and using the Feyman gauge expression 

v liV [x-y) = 5 liV J ^-j , (3.1) 

with 

d 4 p 



* 2tt 



and 



^=4^ sin 2 Qfv) 



we obtain the result that a single loop excitation is suppressed by the factor k exp(— g /4) 
[23]. We have folded in the other factors from (2.27): -1 for a fermion loop, -1 from the 
signed KS phase, and k from the hopping parameters. Now, consider the effect of such a 
loop on a photon propagating between electric currents on the direct lattice at sites and 
y. From the disentangling transformation (2.23) we know that electric current interacts 
with the four-dimensional curl of the Dirac sheet, which for the simple dual plaquette 
excitation is a plaquette on the direct lattice. We can then consider the process shown in 
figure 2. 

Summing over all possible positions and orientations of the loop, we find for the 
contribution v^l to the "dressed" propagator: 

u{S(0,y) = 2 9 2 ^2^2 v /J , a (0,x)x[2v au (x,y) - v av {x + [3, y) - v au (x - $,y) -vp u {x,y) 

x a>/3 

+ vp v (x + a,y) + vp v (x - $,y) - v p „(x + a-$,y)]. 

(3.2) 

It is most convenient to evaluate (3.2) in Landau gauge; 

f exp(ip.(x -y))^ ( , ,„ „v 

Vp V (x,y) = J ^ V„„(p), (3.3) 
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with the lattice transverse projection operator V^ u (p) given by 



^(p)=v- ^HT 1 (i r ) - ^ 



We find 



/• p-ip-y n _ P -iPcx\(-\ _ pip/a'i I 

«£3(o, v) = y 2 y p | ^(p) -E^(p)^(p) ^ I • ( 3 - 5 ) 

In the long wavelength limit p — > the second term in braces is 0(p 2 ) and hence irrelevant. 

Combining this result with the suppression factor of figure 1, we obtain for the full photon 

[o] , [i] 

propagator v^ v = v^l + v^l: 



,-ip.y f f „2 

,2 4 „„„ / y 



lim v^(0,y) = J e -A r V ilv {p) \\ + y ^ 4 exp ^-^jj . (3.6) 

Expression (3.6) is of the same form as the bare continuum propagator, except that 
the overall strength of the interaction has been rescaled, that is, electric charge has been 
renormalised: e R = Z g e 2 , with the renormalisation constant Z g given by 

z » = 1 + 6 exp (-T)' < 3 ' 7 > 

This expression is the first term of an expansion in but is non-perturbative (and 
exact) in g (it is also worth noting that since the fermion determinant is polynomial in 
[i, the expansion in is finite, and there are no non-perturbative contributions). The 
most interesting aspect is the + sign of the correction, which is in contrast to the charge 
renormalisation due to an electric current loop found originally for NCQED in [23]: 

e 2 / e 2 



Z e = l-—^[--y (3.8) 

Of course, the difference has its origin in the factor ieg in the interaction between electric 
and magnetic currents, and supports the claim made in the introduction that monopoles 
anti-screen electric charge. Indeed, at the self-dual point e = g, Z e = Z~ x and bare 
charges remain unrenormalised to this order. This agrees with the results of Calucci and 
Jengo [7], who considered an expansion not in inverse mass but in number of fluctuating 
loops. Because they considered scalar matter fields, they were able to use a proper time 
formalism and consider the current loops as continuum worldlines. They also found that 
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a charge renormalisation was caused by the effect of small intermediate loops (ie. from 
the limit T — > 0, where T parametrises the length of the worldline loop), although the 
precise details of their regularisation are left hazy. It would be interesting to pursue the 
inverse mass expansion further to see if the self-duality persists (indeed, the expansion 
was originally suggested as a means of treating fermion fields in [7]). For larger magnetic 
loop excitations, the dual loop interacting with the electric current is no longer the same 
shape; indeed, a large planar monopole loop corresponds to an array of elementary electric 
current plaquettes in the orthogonal directions. It seems highly plausible, however, that 
Calucci and Jengo's geometrical arguments that electric and magnetic loops of the same 
size and shape will induce opposite charge renormalisations will continue to hold in the 
limit that loop size is much less than photon wavelength, although the fluctuating signs 
due to the fermionic nature of the currents may present a complication. 

Now, in order to define a continuum field theory, we need to approach the continuum 
limit, in which all particles are light; the 1/fx expansion cannot help in this regard. Next 
we consider the limit e — > 0, g — > oo. By self-duality, our findings will also hold in the 
dual limit e — > oo, g — > 0. For g sufficiently large (see below), the exponential term 

2 

exp(— ^-m^v^rn^) in (2.27) alone will suppress monopole polymers, and the monopole 
sector will be described entirely in terms of monomers and dimers, which do not interact 
with the electromagnetic field. In this limit the argument of the third exponential term in 
(2.27) describing the monopole-electron interaction also vanishes; only loops which subtend 
a non-zero solid angle feel this term. Therefore Zqemd factorises: 

lim Zq E md = Zncqed(^)Zncqmd(9)- (3.9) 

e— >0,g— >oo 

By reversing the sequence of transformations leading to eqn. (2.11) we see that the elec- 
tron sector is now governed solely by NCQED (2.1), which may be adequately treated 
in perturbation theory. The staggered lattice fermion formulation used here is known to 
describe four physical fermion species; in the long wavelength limit we recover continuum 
QED with Nf = 4. The monopole sector, on the other hand, is saturated by monomers 
and dimers. It is known rigorously [24] that in the chiral limit fx — > (k — > oo) the chiral 
symmetry of the monopole degrees of freedom is spontaneously broken. In continuum field 
theory, this is signalled by a non- vanishing condensate (ijrij)), where t/> are monopole field 
operators. In the polymer language this symmetry-breaking order parameter is given by 
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[21,24] 

lim (#(,))4 91n ^ MD =k* (3.10) 

Symmetry breaking is signalled by the equilibrium concentration of monomers vanishing 
with some power of p less than or equal to one. 

In the staggered fermion formulation, global chiral transformations are defined by a 
U(l) rotation: 

ip(x) h- > exp(iae(x))ip(x) ; ${x) \— > exp(iae(x))'ijj(x). (3-H) 

Hence in the polymer representation these rotations leave both polymers and dimers 
uneffected, since each contain an equal number of odd and even sites. Electron and 
monopole fields can be independently rotated - so in the p — > limit Zqemd has a 
U(l)y <S> XJ(1) A <S> U(l)y <S> U(l)^ symmetry, with V and A standing for vector and axial 
vector symmetries respectively on both direct and dual lattices. Thus for e — > g — > oo, 
the condensate (N^) breaks the symmetry to U(l) y <8> U(1) A <8> U(l)y. Of course, by 
Goldstone's theorem there must in this case be a massless boson in the spectrum - this 
will be a tightly bound monopole - anti-monopole state which we shall refer to as a "mag- 
netopion". In the g — > oo, \i — > limit in which the monopole physics is described by 
dimer configurations, the magnetopions have vanishing spatial extent, and do not couple 
either to photons or electrons. 

For g large but finite, one may once again estimate the effect of small monopole loop 
excitations on charge renormalisation, using a mean field argument. If the vacuum consists 
of a fraction p of monomers (p = 1 in the 1/p expansion), then the excitation cost of a 
monopole loop changes to (2p) 4p . Using (3.10), the expression for Z g becomes 

Z g = 1 + 2g 2 exp (-^ + 4(^)a» m ^ • (3-12) 

Note that the correction involving (ipip) vanishes in the chiral limit. For electric current 
loops we can use the perturbative QED result 

Z e = l-^ln-. (3.13) 

D7T p 

Since the loop suppression is now due to the factor exp(— g 2 eqn. (3.12) is no longer 
the first term in a systematic expansion; the next smallest loop is suppressed by a factor 
exp(-0.4311^ 2 ) [23]. 
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Before turning to the behaviour of QEMD for intermediate values of e, g, it will 
be useful to discuss and contrast the phase structure of the two other models we have 
introduced, NCQED (2.1), and Villain QED (2.12). As mentioned above, in the long 
wavelength limit NCQED resembles continuum QED with Nf = 4. Numerical simulations 
of this model reveal that in the limit \x — > it exhibits a continuous chiral symmetry 
breaking phase transition for e in the range 2.2 - 2.3 [13,14]. The Villain model, which has 
no fermion fields, has a phase transition, possibly weakly first order, at e ~ 1.25 [25], or 
g ~ 5.0, separating confinement and Coulomb phases. There is no local order parameter 
due to Elitzur's theorem. 

At first sight, it is not clear in this language why either model has a phase transition. 
For large x, v(x) oc 1/x 2 , so that any change in the coupling strength in (2.9) or (2.18) 
can simply be absorbed by a change of length scale. Thus for weak coupling, Villain QED 
resembles continuum QED, and electric charges interact via Coulomb's law. However, 
for small x, the lattice cutoff breaks this scaling symmetry; eg. v(0) is actually a finite 
quantity. This means that in the Villain model, once g is sufficiently small the monopole 
excitations of like sign no longer have a strong repulsive interaction, and those of unlike 
sign are no longer strongly attracted. At this point hitherto small monopole loops become 
free to grow and spread over large distances, forming a plasma which screens the long- 
range Coulomb forces. If a Wilson loop of electric current is introduced into the plasma, 
the interaction term exp(— i2nm^v^, (eAS)^) has the effect of disordering its phase, and 
causing its expectation value to decay as its area, signalling confinement. This is a non- 
perturbative generalisation of the anti-screening effect discussed above. For a loop of 
length L the activation energy ~ g 2 Lv(0) : and the entropy ~ Lin 7 - so the critical g can 
be estimated by an energy /entropy argument [19]. This "monopole condensation" is now 
widely accepted as being the mechanism behind the phase transition in not only the Villain 
model but also U(l) lattice gauge theory with the Wilson action - there are similarities 
with the Kosterlitz-Thouless-Berezinsky mechanism in the two dimensional X-Y model. 

Now consider NCQED (2.9). This model has the same driving term describing a 
Coulomb gas of current loops, but with three important modifications. First, the loops 
are constrained to be singly-charged and self-avoiding. Secondly, there is the sign factor. 
Thirdly, even for a given loop configuration {T} there will be many different background 
monomer/dimer configurations, each carrying its own weight K~ Nm . In the limit e — > oo all 
loops are suppressed, and we are left with a pure monomer/dimer system which is known 
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to break chiral symmetry [24]. In the limit e->0we are left with free massless fermions 
- hence the correlation length must diverge. In this limit large loops are unsuppressed, 
although it is far from clear whether they must dominate Z^cqed - indeed since the 
contributions to Z contain fluctuating signs an entropy argument cannot be constructed. 
Thus, it is unclear whether the phase transition in this model coincides with loops growing 
without bound - so much is hidden in the sign factor. Even continuum QED can potentially 
exhibit singular behaviour, via the unbounded o^ v F^ v term which appears in the expression 
for the inverse propagator [26]. 

Now we return to QEMD; for intermediate values of e, g, the interaction between 
electrons and monopoles governed by the final exponential in (2.27) becomes crucial. As 
discussed above, the standard lore is that monopole loops disrupt the phases of electron 
current loops, causing contributions from larger electric loops to cancel in the partition 
function, and hence suppressing these loops, leading to electric charge confinement. Simi- 
larly, electron loops will disrupt magnetic current loop phases. We may then question the 
existence of a phase in which large polymers are present simultaneously on both direct and 
dual lattices (which would correspond to the coexistence of light electrons and monopoles) . 
Even though the loops each carry a sign from the polymer expansion of the determinant, 
making each term in the partition function of indefinite sign, the additional random phase 
due to the interaction must surely effectively cancel all contributions to Zqemd of this 
form. This physical interpretation is based on the idea that the polymers are the nearest 
we can get to a "fermion worldline". We should also not forget the technical obstacle, 
discussed in section 2, namely a demonstration that Zqemd is positive for real values of 
e. 

To summarise our conclusions so far: there is probably no extended region of the phase 
diagram where both electrons and magnetic monopoles are light propagating particles. For 
e — > the only light degrees of freedom are electrons and decoupled magnetopions, and 
vice versa for g — > 0. For g large but finite, the monopoles' chiral condensate (ipip) is non- 
vanishing in lattice units, and hence divergent in physical units. By analogy with hadronic 
physics, therefore, we expect the decay constant /„■ to diverge and the pion screening 
length m^l ' f% to vanish in the chiral limit, leaving a theory of non-interacting bosons in 
the monopole sector. The question of the magnetopions' coupling to photons and hence 
electrons is more delicate, depending on the spatial extent of the polymer configurations 
dominating Zqemd, but from the previous discussion there is no particular reason to 
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suppose that the coupling will remain non-vanishing in the continuum limit if {(pip) > 
The only scenario which supports an interacting continuum limit for QEMD, therefore, 
is one in which the chiral symmetry of the electrons is spontaneously broken precisely 
at the self-dual point e = g = V2tt; by duality the monopoles' chiral symmetry will be 
restored at precisely the same point. It may thus be possible to approach the self-dual 
point from either phase to take a continuum limit, in which the electrons are in their 
chirally symmetric phase (say), with physical mass vanishing in the limit \i — > 0, and 
the monopoles are in a chirally broken state, with the condensate (ipip) vanishing non- 
analytically as e — ■> (\/27r)_, and some related mass scale defining an inverse correlation 
length, giving dimensional transmutation a la Miranskii [11]. Of course, the self-dual point 
is also a candidate for a zero of the /3-function using the arguments following (3.8). 

Unfortunately, the attractive scenario of a fixed point at the self-dual point appears 
to be excluded by current numerical data. We have already noted that in the limit g — > 
oo the monopole loops decouple to leave standard lattice NCQED. As discussed above, 
NCQED is known to exhibit a chiral symmetry breaking phase transition at e ~ 2.2 — 2.3, 
which is some way short of the self-dual value e = 2.5066 ... It is difficult to believe (and 
indeed runs contrary to the arguments presented above) that the effect of monopole - 
electron interactions will revise this value upwards: if anything monopoles make the e + e~ 
interaction stronger, and would thus reduce the value of the critical coupling needed for 
chiral symmetry breakdown. Hence we must conclude that at the self dual point electrons 
and monopoles probably both have their chiral symmetry spontaneously broken. The 
resulting model is necessarily massive, the correlation length thus finite, and the full j3- 
function no n- vanishing at this point. Instead there must be a region of non-zero width 
centred on the self-dual point, separating two phases where light particles can exist at two 
critical couplings e c , 2iv/e c . 

One final point should be made. As mentioned in the introduction, the considerations 
of this paper were originally motivated by the observation that monopole-like excitations, 
defined on the dual lattice, become dense in a geometric sense and percolate very close 
to or actually at the point where chiral symmetry is broken [13,18,27]. This observation 
encouraged us to speculate that monopole condensation was the agent behind chiral sym- 
metry breaking in NCQED [26]. The truth cannot be so simple - as pointed out in [28], the 
percolation transition does not correspond to a condensation in the Bose-Einstein sense 
(assuming bosonic monopoles), and indeed the monopole excitations of NCQED do not 
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have the long range potentials that could have the effects discussed here. Only formula- 
tions in which the action is periodic in the field strength, and hence in which the Dirac 
string is invisible, can give rise to the monopole condensation phenomenon usually dis- 
cussed. Moreover the monopole-like excitations of [18] will still be present as dislocations 
in the photon field of the NCQED sector of the lattice QEMD model, independent of 
the dynamical monopoles we have introduced explicitly. Therefore it seems unlikely that 
the dynamical monopoles described here could be the primary agent of chiral symmetry 
breaking in QEMD, at least for Nf = 4. It is interesting to note that the numerical co- 
incidence of chiral and monopole percolation transitions at higher values of Nf remains 
impressive [29]. Moreover, the critical exponent v derived from the size of the largest 
cluster is consistent with the critical indices found in powerlaw fits of the chiral equation 
of state of NCQED [29]. Therefore it may well be that the monopoles seen in NCQED 
are connected in some way with chiral symmetry breaking. Quenched simulations using 
alternative actions for the lattice photon field manage to separate the transitions, but still 
suggest that chiral symmetry breakdown is due to fluctuations on the scale of the lattice 
spacing [16]. 

4. Discussion 

We have succeeded in formulating QEMD using an explicit regularisation which pre- 
serves an electromagnetic duality symmetry. This has been done in effect by introducing 
an interaction term into the expansion of the partition function, rather than directly into 
the action. Our formulation (2.27) shares the non-local nature of previous approaches, but 
is more economical in the sense that the Dirac quantisation condition necessary to ensure 
covariance of Green functions is minimal, viz: 

eg = 2n. (4.1) 

This follows from the natural appearance of monopole-like excitations in abelian lattice 
gauge theory [19]. The lattice formulation admits non-perturbative approaches such as the 
1/n expansion discussed in section 3. Although we have not pursued this in detail, it seems 
plausible that higher order calculations will prove to be consistent with the claim of Calucci 
and Jengo that monopole and electron virtual loops give radiative corrections to the photon 
propagator of opposite sign, and hence that there will be no charge renormalisation at the 
self-dual point e = g. 
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On the negative side, the fermionic nature of the monopoles, which we impose pre- 
cisely in order to have a self-dual symmetry, means that it remains to be proved whether 
Zqem d(g) > for e real. Moreover, the fluctuating signs in the terms of Zqemd 
mean that Monte Carlo simulations are probably impracticable; numerical studies of low- 
dimensional fermionic theories using the polymer representation, and even free fermions, 
have proved difficult [21,30]. Despite this, one can argue that results from simulations of 
NCQED suggest that QEMD in the form (2.27) does not have a continuum limit at the 
self-dual point, and that there is no value of e for which both electrons and monopoles are 
light particles. 

The reason for the difficulty of finding a continuum limit is probably the excitation 
energies for polymer loops discussed in section 3. Even the smallest loop, a single plaquette, 
is suppressed by a factor exp(— e 2 /4) ~ 0.21 at the self-dual point - large loops necessary 
for the description of light particles are even further suppressed. The suppression factors 
are intimately bound up in the geometry of the lattice regularisation; it is difficult to see 
how the lattice could be modified without also destroying the duality symmetry. 

There are other possible modifications. Recall that the reason a continuum limit 
appears elusive is that the critical coupling e c of NCQED is somewhat smaller than \/2n. 
Could NCQED be modified to yield a larger e c ? One could introduce a repulsive four-fermi 
interaction to delay the onset of spontaneous chiral symmetry breaking. The suggested 
lattice interaction is [31] 

(4.2) 

A 4 

where x, x are generic fermion fields and 4>^ is a vector auxiliary field, defined on the lattice 
links, which is freely integrated over. The + sign in the definition of M is responsible for 
the repulsion; it renders M hermitian rather than antihermitian, and hence forces the 
eigenvalues of the full fermion kinetic operator to be complex. Repulsive fermion - anti- 
fermion interactions are not usually considered in isolation due to Dyson's argument about 
vacuum stability [32] - in the context of a strongly-coupled theory these arguments are 
hopefully inappropriate. 

A second possibility to increase e c is to increase the number Nf of fermion flavors in 
the model. Simulations with varying Nf [29,33] show clear evidence that e c increases with 
Nf. For sufficiently large Nf (estimated between 13 [33] and 24 [29]) the phase transition 
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becomes first order. However, it is conceivable that there is a window of 8 < Nf < 15 
where both the transition is continuous and e c > \f 7 2m (ie. 1/ej? = 0% = 0.159 . . .). In this 
case introduction of a dual sector of magnetic monopoles could force the existence of a 
phase transition, and hence a fixed point, at the self-dual point. 

Of course, even if no interacting continuum limit can be found, the model presented 
here may be of interest as an effective theory of strongly interacting and massive fermionic 
abelian monopoles arising from some more fundamental model. It is possible in this case 
to relax the condition that ^monopole = ^electron- It is both interesting and amusing 
to contemplate the phenomeno logical implications of lattice QEMD. A recent survey of 
experimental and observational limits on monopole matter, and in particular the effect 
of virtual monopole loops, has been given in [34]. The novelty of the present approach 
is that monopoles are strongly interacting fermions, which means that in our world of 
small e ~ 0.3, monopole chiral symmetry is spontaneously broken, and there will be at 
least one light Goldstone particle in the monopole sector, the magnetopion ix m . Can 
magnetopions interact with ordinary matter? The scattering of photons off tiny magnetic 
dipole moments will depend on the details of the magnetic charge distribution in the 
magnetopion, which in turn depends on the (non-perturbative) dynamics. However, we 
note that magnetopion production via virtual photon decay 7 — > mv m is forbidden by 
charge-conjugation symmetry, just as strong interaction decays such as p° — > utt are 
forbidden. The main interaction of magnetopions, just as for axions, will be via the 
coupling 71^77; indeed, massive magnetopions will decay to two photons. Since the lattice 
regularisation presented here does not admit an axial anomaly, magnetopion decay will be 
suppressed by a factor of 0(m^ ) by the usual arguments [35]. This raises the interesting 
question of whether there is a range of bare mass \i such that the magnetopion is both 
sufficiently long-lived and sufficiently heavy to be of cosmological significance. Such issues 
cannot be resolved without a further understanding of the detail of the strongly interacting 
monopole dynamics. 
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Figure 1 

Photon corrections to the smallest loop excitation 




0; jLt 

Figure 2 

Loop correction to the photon propagator 
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